Exchange of stability as a function of system size in a nonequilibrium system 
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In equilibrium systems with short-ranged interactions, the relative stability of different thermo- 
dynamic states generally does not depend on system size (as long as this size is larger than the 
interaction range). Here, we use a large deviations approach to show that, in contrast, different 
states can exchange stability as system size is varied in a driven, bistable reaction-diffusion system. 
This striking effect is related to a shift from a spatially uniform to a nonuniform transition state and 
should generically be possible in a wide range of nonequilibrium physical and biological systems. 



The statistical physics of nonequilibrium systems has 
proven to be an enduring source of unexpected and in- 
triguing phenomena. Historically, studies of driven sys- 
tems have tended to focus on the thermodynamic limit of 
infinite system size. In recent years, however, experimen- 
tal advances have made it possible to study everything 
from micromagnets and nanotubes to living cells at the 
mesoscopic scale, and this in turn has led to a growing 
theoretical interest in finite-sized stochastic systems. A 
major problem is to understand the relative stability of, 
and transition rates among, different (mcta)stablc states; 
this issue has been studied extensively for some equilib- 
rium, detailed-balancc-obeying models [U-H)]. Here, we 
consider the same question in a simple model of a far- 
from-equilibrium, driven chemical system and show that 
far richer behavior is possible when detailed balance is 
violated. In particular, we find that — unlike in an equi- 
librium system — the system's two stable states can ex- 
change stability as the system size increases. This rever- 
sal is related to competition between a homogeneous and 
a spatially- varying transition state [HIH, an d we expect it 
generically to be possible whenever there is no underlying 
Boltzmann distribution to ensure that the relative stabil- 
ity is independent of the transition path. Similar results 
may thus apply to a wide range of systems of interest in 
condensed matter, chemical, and biological physics. 

In what follows, wc start by introducing the Schlogl 
model of a bistable chemical system and formulating its 
mean-field rate equations in terms of a deterministic po- 
tential V(c). The minima of this potential correspond 
to two locally stable states, and one might naively imag- 
ine that the state with the lower V is the more stable 
of the two. Indeed, in a spatially extended system with 
diffusive transport, the mean-field equations have front 
solutions in which the state with lower V always invades 
that with higher V. In a well-stirred system, in contrast, 
only noise-induced transitions between states are possi- 
ble, and the mean-field description provides no informa- 
tion about relative stability. Surprisingly, in this case, an 
analysis that accounts for the fundamentally stochastic 
nature of chemical reactions shows that, in certain pa- 
rameter ranges, the state with higher V is found with 
greater probability 0, 0] • To reconcile these two seem- 
ingly contradictory predictions, we consider a model in- 



cluding both diffusion and fluctuations. We first show 
numerically that the two states exchange stability as the 
system's spatial size grows. We then use a semiclassical 
approach to explain this dependence. By placing bounds 
on the action that determines the transition rate between 
states, we are able to show that the relative stability 
must follow the deterministic potential for large enough 
systems, where the transition occurs through nucleation 
followed by deterministic front motion, while the well- 
stirred result applies for small enough systems. 

The bistable Schlogl model [t| consists of the chem- 
ical reactions 
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2X + B 



(1) 

where the concentrations of A and B are held constant. 
The corresponding mean-field rate equation for the con- 
centration c of X in a well-stirred system is 

c = (fco + he 2 ) - (k 3 c 3 + he) = -V'(c) , (2) 

where we have absorbed the concentrations of A and B 
into the rates h and h, respectively, and the prime de- 
notes a derivative. This equation defines the determinis- 
tic potential V(c). If we move away from the well-stirred 
limit and let c depend on a spatial coordinate z, a diffu- 
sion term must be added, and Eq. @ generalizes to 



dic= -V'(c) + Vd 2 z c 



(3) 



For appropriate choices of the h, V(c) has one local 
maximum c s and two local minima c\ and C2, with c\ < 
C2, corresponding to two (meta)stable states. Eq. |(2J) is 
invariant with respect to a simultaneous rescaling of con- 
centrations and of rate constants, a fact we can exhibit 
explicitly by introducing a typical concentration scale Cq 
of the same order as c\ and c% and writing the h as 
h = Ai(co) 1-1 , where the A, have dimensions of inverse 
time. For fixed Xi, the dynamics of x — c/cq is indepen- 
dent of Co. 

The existence of two (meta)stable states invites the 
question of their relative stability. One expects that, 
when noise is properly taken into account, c will be found 
with high probability near c\ or C2, and one might guess 
that V(c) determines which of the two is more probable 
(i.e. stable) and which less probable. Indeed, in the limit 



2 



of infinite system size, Eq. ([3j admits traveling front so- 
lutions of the form c{z, t) = f(z — vt) in which the state 
with lower V expands into the one with higher V 0, 0, HI ■ 
Eq. ([3]), however, is only a mean- field approximation 
to a more realistic model that accounts for the random 
nature of the molecular collisions that lead to chemical 
reactions. To incorporate both these intrinsic stochastic 
effects and diffusion, we use a mesoscopic compartment 
model. In general, a d-dimensional reaction vessel is par- 
titioned into M elementary compartments of linear size 
h and volume V = h d , and each molecule can jump be- 
tween neighboring compartments but can react only with 
other molecules in the same compartment. The station- 
ary probability P(n 1 , n 2 , . . . , n l , . . .) to find ri 1 molecules 
in compartment i then satisfies 

[W+(n* - 1)P(., n l - 1, .) - W+{rf)P{., n\ .)] + 

i 

Y [W-(n l + 1)P(., rj + 1, .) - W-(n l )P(., n\ .)] + 

i 

D Y [(n 1 + 1)P(., n* + 1, n' - 1, .) - n*P(., n<, n', .)] 

-0, (4) 

where D is the jump rate between compartments and the 
last sum is taken over neighboring compartments 
The reaction rates are 

W+(n) = k V+(k 2 /V)n(n-l) , 

W_(n) = k 1 n+(k 3 /V 2 )n(n-l)(n-2) . (5) 

In general, this model allows for barriers that slow diffu- 
sion between compartments, but we are primarily inter- 
ested in using it as an approximate description of a spa- 
tially continuous system. One can show that the average 
behavior of the stochastic compartment model (j4]) ap- 
proaches that of the continuum reaction-diffusion equa- 
tion ©, with V = h 2 D, when Dt r = Vt r /h 2 > I and 
cqV 3> 1 [l(J. Here t r ~ 1/Aj is a typical time be- 
tween reactions of an individual molecule. The first in- 
equality ensures that individual compartments are well 
mixed; while this could be accomplished by decreasing h 
at fixed T>, the second inequality demands that this not 
be done at the expense of having very few molecules in 
each compartment. 

Although there is no general analytical expression for 
the stationary distribution satisfying Eq. (j4|) , one can be 
found when M = 1 and the entire reaction volume is well 
stirred. In this case the total number of molecules n is 
distributed according to 



Pws(n) 



w+(j-i) 

W-(j) 



(6) 




FIG. 1: Logarithms of the probability distribution P wa (x) 
in the well stirred model and of the deterministic potential 
V(cox) Both quantities are scaled by their values at xi. 



that the typical number of particles per compartment 
Q, = CqV becomes large, and it is thus useful to rewrite 
P ws in terms of x = c/cq = n/Q as 11 1 



P ws (x) = K(x;Sl)e- as ™W 



(7) 



where K(xi\Q,)/K(x2',£l) is bounded as £1 — > oo for ev- 
ery xi, X2 > 0, so that exp[— fi5 ws (x)] determines the 
dominant large £1 contribution to any probability ratio. 
The action S ws (x) is defined by 



dS ws (x)/dx = In [w- (x) / w + (x)] 



(8) 



where K is a normalization constant. Below, we will be 
particularly interested in asymptotic results in the limit 



where w + {x) = Ao + X-ix 2 and u>_ (x) = X\x + X3X . It is 
clear from this equation that the extrema of S ws (x) and of 
V(cqx) occur at the same values of x. The two functions, 
however, can otherwise be very different. Indeed, for the 
parameters used in Fig. Q] [l2j], V(c2) — V(ci) and A ws = 
S W s(x2) — S ws {xi) (where the extrema Xi = Cj/co) have 
opposite signs. Since Eq. J7]) implies that it is the sign 
of A ws that determines which of the two states is more 
probable, the deterministic potential does not reliably 
predict relative stability in the well-stirred case. 

Semiclassical approach. In a well-stirred system, 
no transitions between states can occur without fluctu- 
ations. In contrast, the kinetics ([3]) with diffusion does 
allow for front-driven transitions between states even in 
the mean-field limit. This suggests that, unlike in the 
well-stirred system, there may be cases in spatially ex- 
tended systems where the deterministic potential V(c) 
does in fact determine relative stability. Direct simula- 
tions of the fully stochastic compartment model (using a 
kinetic Monte Carlo algorithm with separate treatment 
of reaction and diffusion steps Q) demonstrate that this 
intuition is correct. Indeed, Fig. [2] shows that the two 
states can exchange stability as M is increased; for small 
M, their relative stability is the same as in the well stirred 
limit, but for larger M, the state with lower V regains the 
upper hand. As we now show, this stability inversion can 
be understood within a semiclassical approximation. 

Just as in the well-stirred case, the stationary distri- 
bution P satisfying Eq. Q can be written for large f2 
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FIG. 2: The distribution of the average scaled concentration 
x = £>j/(Afn) for different M and V = 20 [H. /nset: 
schematic of the compartment model, with molecules allowed 
to hop between compartments and to react within each com- 
partment. 
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P{x\ 



) = K M (x\. 



(9) 

where x 1 = and D, = cqV remains the typical number 

For well-mixed sys- 
is often thought of as 
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of particles per compartment 
tems, the eikonal approximation I 
a large volume approximation, but here the requirement 
that the compartment model approximate a continuous 
reaction-diffusion system constrains the volume of each 
compartment, and it is more useful to think of making f2 
large by lettin g en — > oo with the Aj, D, and V fixed. 



It is known [llj, [13J that the function S(x) can be ex- 



pressed as the minimal action 

S(x) — S(xi) = minmax < / dt Ipx — H (p, x)] 

H(p,x) 

J2(e pt - 1) \w+(x l ) - w-(x i )c- pi +DA l (xc- p ) ,(10) 



where stable stationary point of the deterministic 

kinetics, x a point in its basin of attraction, and Aj the 
appropriate discrete Laplacian centered on compartment 
i. The extrema of S(x) correspond to the fixed points of 
the dynamical system 



x 1 = w + (x l ) 



W- 



\x l ) + DA t (x) 



(11) 



and have the same stability [ll| . 

In the continuum limit h — > 0, Eq. ([lip becomes a 
rescaled version of Eq (|3l), whose stationary states are 
well-characterized Q, Il4ll6j . It is reasonable to ex- 
pect (see Supplemental Material [SM]) that the station- 
ary states of the discretizcd version have similar prop- 
erties for small enough h. Then, with reflecting or 
periodic boundary conditions, only the uniform states 



Xi = (x\, . . . , x\) and X2 = (x2, ■ ■ ■ , x%) are stable, and 
there is a unique (up to symmetries) fixed point x s of (fTTj) 
with only one unstable direction. This saddle is used to 
define the function S(x) everywhere, i.e. to fix the value 
of S(xi) and thus of the stability index A through a 
matching procedure [III, EH : 

A = S{x 2 ) - S( Xl ) = AS hs - AS 2 , S , 



AS\ s = minmax / dt Ipx — H{p, x)\ , 

a»(*) P(*) J Xl 
[ x » 

AS-? s = minmax / dt \px — H(p, x)] . 
»(*) Pit) J X2 



(12) 



11] 



Both A5j iS are non-negative 

We now focus on the regime where (as in Fig. [T]) 
V(c2) — V(ci) < 0, but A ws > 0, and study the sign 
of A, and through it the relative stability of the two uni- 
form states. This sign strongly depends on the saddle x s , 
which, in the continuum limit — >■ 0, is spatially uniform 
for a small enough system but becomes nonuniform at a 
critical linear system size [l], 0, HI • In the limit of large 
system size, the saddle profile is close to ii everywhere 
except in a localized region whose size remains constant 
as the system size grows. This form reflects that fact 
that the deterministic traveling waves favor the state x 2 , 
so that only a small nucleus is required to initiate a tran- 
sition from uniform xi to uniform x<i- Similarly, for our 
compartment model, the stationary states and optimal 
trajectories in Eq. (fT2"]l arc uniform for M less than some 
M c , and one can show that A = A/A ws . The stability is 
then that of a well stirred system, independent of D. For 
M > M c , solving the double optimization problem (fl"2"|) 
is extremely difficult. One can, however, verify numeri- 
cally (SM, Fig. @| that x s has the same shape as in the 
continuum limit, and thus deviates appreciably from X\ 
only in a small region. Using this fact, one can derive 
bounds on A. 

In particular, it is not hard to see that A5i iS , which 
describes the difficulty of reaching x s from x\, remains 
smaller than a fixed constant, independent of M. This 
bound is simply given by a particular trajectory solving 
the first maximization problem in Eqs. (|12p . As x Sl and 
thus the trajectory, differ from x\ only in a small region, 
the total action is finite, even for infinite M (see SM). 

Similar reasoning indicates that AS 2 . S should grow lin- 
early with M. Indeed, in the limit M — > oo, the region 
where x s differs from X\ becomes negligible, and one can 
focus on an optimal trajectory that must take the system 
from x 2 to x\ on an essentially infinite domain. Physi- 
cally, one expects that such an optimal trajectory should 
correspond to a front solution of the variational equations 
traveling with constant speed. Each successive compart- 
ment then makes the same contribution to the action, 
and one has AS 2 , S (M) ~ MAq for some Ao- (See also 
SM.) Moreover, it is clear that for any integers M and I, 
AS 2 , S (M) > AS 2 ,s(Ml)/l, because an optimal trajectory 
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FIG. 3: The scaled action —S/M for different numbers of 
compartments. The plotted values represent the logarithm of 
the probability of an average concentration x — '^2 i x % /M in 
Monte Carlo simulations [Sj], extrapolated to Q, = oo. Curves 
are shifted such that S{x\) = in order to clearly indicate the 
stability exchange. Note that the peaks in the curves move 
towards x\ as M grows, reflecting the fact that most of the 
saddle profile is near x\ for large M. 



for a system of size M can be mirrored I times to create a 
suboptimal trajectory for a system of size Ml. Thus, the 
asymptotic large M behavior AS , 2 jS (M) ~ MAq implies 
that AS 2 , S {M) > MA even for finite M. Since AS 1>a 
remains finite as M increases while AS^s grows without 
bound, A = AS\ tS — AS^s must change sign at some 
71/ > M c , and the two states exchange stability as the 
system size increases (Fig. 

Discussion. The large deviations approach that we 
have used to show that an exchange of stability must 
occur is formally valid in the limit that f2 —y oo at fixed 
compartment volume V and number of compartments M. 
This order of limits is important. In particular, because 
the prefactor Km in Eq. depends on M, it can over- 
whelm the exponential factor if M tends toward infinity 
at fixed f2; thus, our calculations describe transitions be- 
tween (meta)stable states in a finite-sized, mesoscopic 
system but cannot be used to study the phase transition 
in the Schlogl model defined in the limit of infinite sys- 
tem size [17] . Similarly, our results hold in the limit of 
a large number of reacting particles and thus describe 
a distinct phenomenon from the exchange in the most 
probable state observed in a well stirred system with a fi- 
nite (small) particle number [lj| or the recently-described 
noise- induced reversal of front propagation direction Q . 
Our calculations also go considerably beyond a previous 
study that examined transition rates very near a bifur- 
cation in the Schlogl model but did not consider relative 
stability [H . Finally, although the exchange of stability 
discussed here could not happen if the transition state 
did not become non-uniform at M = M c (leading to a 
non-analyticity in A) 0,13], the actual exchange occurs 
at some M > M c and need not coincide with any further 
bifurcations of the transition state. 

Although we have presented our results for a specific 
model of a bistable chemical system, the system-sizc- 



dependent relative stability that we describe is far more 
general and should in principle be possible in models de- 
scribing everything from pattern formation (l9j to ecolog- 
ical population dynamics [H or the complex reaction net- 
works present in living cells [20| ■ All that is absolutely re- 
quired is multistability and violation of detailed balance. 
In some cases, more complex stability diagrams are likely 
to be possible. For example, if we no longer insist that 
the compartment model studied here be a good approx- 
imation to a continuum system and thus allow ourselves 
to vary D arbitrarily, we expect a reentrant exchange of 
stability: As D —> oo, the entire system should be well 
mixed, while each individual compartment behaves as a 
separate well mixed vessel as D — > [2l|. In either limit, 
the relative probability of the two states should depend 
on A ws , while for intermediate values of D the sign of A 
may differ from that of A ws . These and similar effects 
should be accessible within the same formalism employed 
here. 

Acknowledgments. We are grateful to Charlie Do- 
ering, Baruch Meerson, and Len Sander for helpful con- 
versations. This work was funded in part by NSF grant 
DMR1056456. 



[1 

[2 

[3, 

[1 
[5 

[6 

[7] 
[8] 

[9 
[10 

[11 
[12 

[13 

[U 

[15 
[16 

[17 
[18 

[19] 



W. G. Faris and G. Jona-Lasinio, J Phys A 15, 3025 
(1982). 

R. S. Maier and D. L. Stein, Phys Rev Lett 87, 270601 
(2001); Proc SPIE 5114, 67 (2003). 

J. Burki, C. A. Stafford, and D. L. Stein, Phys Rev E 
77, 61115 (2008). 

B. Meerson and P. V. Sasorov, Phys Rev E 83, 22 (2011). 
P. Hanggi, H. Grabert, P. Talkner, and H. Thomas, Phys 
Rev A 29, 371 (1984). 

C. R. Doering, K. V. Sargsyan, L. M. Sander, and 

E. Vanden-Eijnden, J Physics-Cond Mat 19 (2007). 

F. Schlogl, Z Phys 253, 147 (1972). 

E. Khain, Y. T. Lin, and L. M. Sander, Europhys Lett 
93 (2011). 

J. Elf and M. Ehrenberg, Systems Biology 1, 230 (2004). 
S. A. Isaacson and C. S. Peskin, SIAM J Sci Comput 28, 
47 (2006). 

M. I. Dykman, E. Mori, J. Ross, and P. M. Hunt, J Chem 
Phys 100, 5735 (1994). 

For all figures, the parameters are ko — .2, k\ — 2, &2 = 
3, kz = 1, and D = 1, in units in which co and A3 have 
been set to unity. 

M. Heymann and E. Vanden-Eijnden, Comm Pure Appl 
Math 61, 1052 (2008). 

M. Buttiker and H. Thomas, Phys Rev A 24, 2635 
(1981). 

E. Scholl, Z Phys B 62, 245 (1986). 

L. Hui, K. Dimitri, and I. Francesco, J Phys A 33, 6413 
(2000). 

P. Grassberger, Z Phys B 47, 365 (1982). 

W. Ebeling and L. Schimansky-Geier, Physica A 98, 587 

(1979). 

U. Bisang and G. Ahlers, Phys Rev Lett 80, 3061 (1998). 



5 



[20] J. Munoz-Garcia and B. N. Kholodenko, Biochem Soc 

Trans 38, 1235 (2010). 
[21] N. Berglund, B. Fernandez, and B. Gentz, Nonlinearity 

20, 2551 (2007); ibid. 20, 2583 (2007). 



6 



SUPPLEMENTAL MATERIAL 

As in the latter part of the main text, we assume here 
that the parameters are such that V(c2) < V(ci) but 
A ws > 0, so that the concentration C2 is favored in the 
mean-field limit while c\ is more stable for a well-stirred 
system. 

The stationary states 

The stationary states of Eq. ^ satisfy 

Vd\c = V(c) (13) 



1.4r 




Compartments 



and represent trajectories of a Hamiltonian system in the 
classical potential — V(c), where the spatial variable z 
plays the role of time. In addition to the constant trajec- 
tories c(z) = c\,c s ,C2, nonuniform trajectories are also 
possible for large enough system lengths. 

These trajectories represent one or more 
(half)oscillations around the local minimum c s of 
— V. Intuitively each swing of these trajectories (or 
every full oscillation for periodic boundary conditions) 
can be independently destabilized, making the number 
of swings equal to the number of unstable directions 
of the corresponding stationary state. Therefore, when 
non-constant stationary states are possible, the saddle 
with only one unstable direction corresponds to a 
trajectory with a half (full) oscillation for reflecting 
(periodic) boundary conditions [2, 14 - 16 1 . 

Fig. [4] shows numerically obtained saddle profiles for 
the discretized model of Eq. (fTTj) . A similar structure 
is obtained as in the continuum limit. In particular, for 
large M the region where the concentration deviates from 
Ci remains the same size even as M increases. (Saddle 
profiles were obtained by solving the coupled algebraic 
equations defining the stationary state by standard meth- 
ods, using the continuum solution as an initial guess, then 
verifying numerically that the linearized dynamics about 
the stationary profiles obtained in this manner had only 
a single positive eigenvalue.) 



Upper bound on ASi tS 
In order to construct an upper bound on the quantity 



ASi s = minmax / dt \px — H(p,x)] , (14) 

*W P(<) J Xl 

it is enough to find a particular trajectory x p (t) between 
Xi and x s and then to evaluate 



FIG. 4: The saddle profile x s representing a fixed point of 
Eq. for periodic boundary conditions and different num- 
bers M of compartments. (The horizontal axis is rescaled 
so that profiles for different M fit on the same graph.) For 
M < M c = 6 the saddle profile is uniform. As the number of 
compartments is increased for M > M c , the profile converges 
to a form where most of the profile, except for a critical nu- 
cleus, is close to the metastable state ci. 



for that trajectory. For a given choice of x p , the maxi- 
mum in this equation is attained for momenta p p satis- 
fying 



= x. 



Because the function 



dH 



(16) 



H(p,x) = 

^(e p! - 1) \w+(x l ) ~ w-(x l )c- p ' + DA x (xc-P)\ (17) 

i 

is convex, Eq. ([Trl)) determines, for each x p , a unique p p 
where a global maximum over all possible p(t) is attained. 

We can take advantage of this uniqueness property 
to specify our particular trajectories x p and p p simul- 
taneously. For future convenience, we rewrite the mo- 
menta in terms of a (so far unknown) function x*(t) as 

log f^fj ■ Knowing that H(p, x) = along the true 
extremal trajectory (llj . we then define our particular 
trajectory by demanding that each term in the sum 

H(x p ,p p ) = 
- 1) [w + (xl) - w. (4)4/4 + DAi(x*)] ,(18) 



P 



S[x p (t)] = max / dt [px p — H(p, x p )j (15) 
pM Jan 



separately vanish: 

= w+(x p ) 



^-^-xl + DAtixJ (19) 
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for each i. These M equations, together with the M 
Hamilton's equations (|16[) . determine x p and p p . In par- 
ticular, we can rewrite (fTB|) explicitly as a set of dynam- 
ical equations for the x l p : 

x p = w + {x p )x p /x\ - w-{x p )x\/x l p + 
D[xJ x p /x* - xjx^ ) + 
D(x^ x p /x t — x p x^/x^ ) , 

where Eq. (JTHJ) implicitly determines the x\ at each time. 
It is not hard to verify that these dynamical equations for 
x p (t) share the same fixed points as the original mean- 
field kinetic equations for the concentrations x(t); indeed, 
for M < M c , where the concentrations remain uniform 
along the optimal trajectory, x p (t) coincides with the 
true optimal trajectory. In particular, one can easily 
check that the uniform state Xp — X 1 ^ — X\ iScL fixed 



point of Eq. (|2QJ for any M. 

We can now use the particular solution (p , x p ) to show 
that AS\, S , remains finite even as M — > oo. First note 
that H(p p , x p ) = by construction. Thus, our bound on 
ASi 8 reduces to 



VpdXp 



A5i, s < / dtp p Xp= / p p dxp = 



(20) 

For concreteness, assume that the single peak in x s is 
centered at i = 0. Away from this peak, x\ — x\ is small, 

and we expect that each integral f** p l p dx l p is similarly 
small. In fact, by linearizing Eq. ([lip about x\ and set- 
ting x l = 0, one immediately concludes that x\ decays 
exponentially to x\ as \i\ —> oo. As long as the maximum 
over t (or cquivalently over x p ) of |pi(t)| approaches a fi- 
nite value in this limit, the sum of integrals with respect 
to the x p is then bounded by a convergent geometric se- 
ries, and ASi^ must remain below a finite bound for all 
M . It is not hard to see that \p p \ in fact cannot diverge 
for large \i\: Eq. (| is linear in cc* and approaches a 
screened diffusion equation for large where x p x\. 
Thus, for large x\ can be decomposed into an expo- 
nentially decaying part whose prefactor depends on non- 
linear behavior in the peak near i = and terms of order 
at most x^. Consequently, neither xl nor p p diverges for 
|i| — > oo, and ASi jS approaches a constant for large M, 
as asserted. 

An interesting remark is that a similar construction for 
AS2, S yields a dynamics analogous to Eq. (|TT1) . but for 
trajectories between x-i and x s , that supports a propagat- 



ing front solution for large enough systems. This estimate 
thus gives an upper bound on AS^.s that grows linearly 
with M, complementing our arguments on a lower bound 
in the main text and in the next section. 

Lower bound on A5*2, a 

In this section, we elaborate on the arguments leading 
to the conclusion AS , 2. s (Af) > MAo for some positive 
Ao, which we obtained in the main text. 

For states x s and X2 connected by a solution x(t) of 
the deterministic kinetics (|11|) . one can rewrite AS2.S as 



AS 2 , S (M) = S(x 2 )-S(x s ) 



»=i 



dS(x) . q 

dx l 

dS(x) 
dx l 



x [w+(x l ) -w-(x*) +DAi(x)] 



(21) 



For large number of compartments M, the solution x{t) 
of the deterministic dynamical system 



x* = w+(x l ) - w-(x l ) + DAi{x) 



(22) 



is a moving front (except in the immediate vicinity of 
the single peak in x s ), and each term in the sum over 
i should thus make the same contribution to the inte- 
gral. One might thus naturally be led to conclude that 
AiS2,s grows linearly with M for large M. One flaw in 
this reasoning is that we have not ruled out the possi- 
bility that dS(x)/dx % depends on M. It will not do so 
as long as the dependence of the action S(x) on the x l 
is sufficiently local, so that the derivative with respect 
to a given x z doesn't "know" how many compartments 
there are in total. Our assumption in the main text that 
the solution to the full variational equations is a mov- 
ing front can thus be replaced by the locality assumption 
that dS(x)/dx l is independent of M for large M. For 
this to be the case, it is sufficient that the second deriva- 
tive gjg^j have only short-ranged dependence on i — j : 

gjg^j < Ae~ c ^~^ for some positive constants A and C 
and \i — j\ large enough. This is known to be true in an 
equilibrium system, but remains to be proven completely 
rigorously for our noncquilibrium model. 



